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Using an isospin-dependent transport model, we study isospin effects on two-nucleon correlation 
functions in heavy-ion collisions induced by neutron-rich nuclei at intermediate energies. We find 
that these correlation functions are sensitive to the density dependence of nuclear symmetry energy, 
but not to the incompressibility of symmetric nuclear matter and the medium dependence of nucleon- 
nucleon cross sections. This sensitivity is mainly due to effects of nuclear symmetry energy on the 
emission times of neutrons and protons as well as their relative emission sequence. We also study 
the variations of the symmetry energy effects on nucleon-nucleon correlations with respect to the 
impact parameter, incident energy, and mass number of heavy ion collisions. 
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I. INTRODUCTION 

Heavy ion collisions with radioactive nuclear beams 
provide a unique opportunity to study the properties of 
nuclear matter at the extreme condition of large isospin 
asymmetry 

H H H H In these collisions, the dy- 
namics is affected not only by the isospin-independent 
part of nuclear equation of state (EOS) but also by the 
isospin-dependent part, i.e., the nuclear symmetry energy 
E sym (p), and the isospin-dependent in- medium nucleon- 
nucleon (JV-JV) cross sections. Knowledge on the den- 
sity dependence of nuclear symmetry energy is essen- 
tial for understanding both the structure of radioactive 
nuclei M, 11 H H| and many key issues in astrophysics 
H E3 Elir Although the root -mean-squared radius of 
neutron distribution in a heavy nucleus and that of a 
neutron star differ by about 18 orders of magnitude, the 
properties of bot h sy stems are strongly affected by the 
same E sym (p) U\m\m Moreover, the nucleosynthesis 
in pre-supernova evolution of massive stars, mechanisms 
of supernova explosions, cooling rates of protoneutron 
stars and associated neutrino fluxes, and the kaon con- 
densation as well as the hadron to quark-gluon plasma 
phase transitions in neutron stars all rely critically on the 
nuclear symmetry energy |9L llCtll4j . However, our knowl- 
edge on E sym (p) is still rather poor despite extensive the- 
oretical studies based on various many-body theories 0] . 
In fact, it is often regarded as the most uncertain prop- 
erty of a neutron- rich matter |15l Il6| . Fortunately, ra- 
dioactive beams, particularly the very energetic ones to 
be available at the planned rare isotope accelerator (RIA) 
and the new accelerator facility at German Heavy Ion 
Accelerator Center (GSI), provide a great opportunity 
to study the E sym (p). It is thus important to investigate 
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theoretically what experimental observables can be used 
to extract information about E sym (p) . In this regard, sig- 
nificant progress has been made recently. For instance, 
measurement of neutron skins of radioactive nuclei via 
their total reaction cross sections [j| or that in stable 
heavy nuclei, such as 208 Pb, via parity-violating electron 
scatterings Q was recently proposed as a sensitive probe 
to E sym (p). In heavy- ion collisions induced by neutron- 
rich nuclei, the pre-equilibrium neutron/proton ratio |17| . 
isospin fractionation 0,[ljJ,|2jj,|2l],|22, isoscaling in mul- 
tifragmentation j23| , proton differential elliptic flow [24| , 
neutron-proton differential transverse flow [25j as well as 
the 7r~ to 7r + ratio 26] have been identified as promising 
probes. 

Observables mentioned in the above involve only the 
momentum, charge, and mass distributions of individual 
particles in the final state of a reaction. On the other 
hand, two-particle correlation functions, through final- 
state interactions and quantum statistics effects, have 
been shown to be a sensitive probe to the temporal and 
spatial information about the reaction dynamics in heavy 
ion collisions at incident energies ranging from interme- 
diate I2H2H2H to RHIC energy 0. In particular, the 
correlation function of two nonidcntical particles from 
heavy ion collisions has been found to depend on their rel- 
ative space-time distributions at freeze out and thus pro- 
vides a useful tool for measuring the emission sequence, 
time delay, and separation between the emission sources 
for different particles 0111 III Hill III . Indeed, time 
delays as short as 1 fm/c and source radius differences of 
a few tenth fm have been resolved experimentally from 
correlation functions of nonidentical particles |3^ |. 

In the present work, we study systematically how 
E sym (p) affects the temporal and spatial structure of 
reaction dynamics in heavy-ion collisions induced by 
neutron-rich nuclei at intermediate energies. A brief re- 
port of present work can be found in Ref. 38]. Wc 
show that average emission times of neutrons and pro- 
tons as well as their relative emission sequence in heavy 
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ion collisions are sensitive to E sym {p). A stiffer density 
dependence of nuclear symmetry energy leads to faster 
and nearly simultaneous emissions of high momentum 
neutrons and protons. Consequently, strengths of the 
correlation functions for nucleon pairs with high total 
momentum, especially neutron-proton pairs with low rel- 
ative momentum, are stronger for a stiffer E syul (p). This 
novel property thus provides another possible tool for ex- 
tracting useful information about E sym (p). In addition, 
we also study how two-nucleon correlation functions de- 
pend on the isospin-independent part of nuclear matter 
EOS, the in-medium JV-JV cross sections as well as the 
impact parameter, incident energy, and masses of collid- 
ing nuclei. 

The paper is organized as follows. In Section [n] 
we discuss the density dependence of nuclear sym- 
metry energy and the EOS of neutron-rich nuclear 
matter. The Isospin-dependent Boltzmann-Uchling- 
Uhlcnbcck (IBUU) transport model used in present study 
is briefly described in Section ITTll In Section llVl we pre- 
sented results from the IBUU model on the nucleon emis- 
sion functions and their isospin dependence in heavy-ion 
collisions induced by neutron-rich nuclei at intermediate 
energies. These include the nucleon emission times, the 
size of nucleon emission source, and the momentum dis- 
tributions of emitted nucleons. In Section ^ we study 
nucleon-nucleon correlation functions and their depen- 
dence on the nuclear symmetry energy. We also present 
results on the variations of the symmetry energy effects 
on nucleon-nucleon correlations with respect to the nu- 
clear isoscalar potential, nucleon-nucleon cross sections, 
as well as the impact parameter, incident energy, and 
masses of heavy ion collisions. Finally, we conclude with 
a summary and outlook in Section IVII 



II. NUCLEAR SYMMETRY ENERGY 



Many theoretical studies (e.g., [3J, |4fl, |4jj, |42j, |43j, |4J ) 

have shown that the EOS of asymmetric nuclear matter 
can be approximately expressed as 



E(p,S)=E(p,S = Q) + E sym (p)6 



(1) 



where p = p n + p p is the baryon density; 5 — (p n — 
Pp)/(Pp + Pn) is the isospin asymmetry; and E(p, 5 = 0) 
is the energy per particle in symmetric nuclear matter. 
The bulk symmetry energy is denoted by the so-called 
symmetry energy coefficient E aym (p) = E(p,S = 1) — 
E(p,S = 0). Its value at normal nuclear matter density 
Po, i.e., E syni (po), is predicted to be about 34 ± 4 MeV 
and is comparable to the value extracted from atomic 
mass data j45j- 

The symmetry energy coefficient E sym (p) can be ex- 
panded around normal nuclear matter density as 
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where L and K syul are the slope and curvature of sym- 
metry energy coefficient at normal nuclear density, i.e., 



L = 3p ( 



dE sym (p) 

' dp lp=po ' 



n ^ 2 G> 2 ff sym (p) , 

9p ° Wp 1 



(3) 
(4) 



The L and K syra characterize the density dependence of 
nuclear symmetry energy around normal nuclear matter 
density, and thus provide important information about 
the properties of nuclear symmetry energy at both high 
and low densities. They are, however, not well deter- 
mined either theoretically or experimentally. The theo- 
retically predicted values for K sym vary from about —700 
MeV to +466 MeV (e.g. [iff), while data from giant 
monopole resonances indicate that K sym is in the range 
from -566 ± 1350 MeV to 34 ± 159 MeV depending on 
the mass region of nuclei and the number of parameters 
used in parameterizing the incompressibility of finite nu- 
clei |47J. There are no experimental data for the value 
of L, and different theoretical models give its value from 
-50 to 200 MeV 0. 

Predictions for the density dependence of E sym (p) are 
thus divergent. In the present study, we adopt the param- 
eterization used in Ref. 48] for studying the properties 
of neutron stars, i.e., 



E synl (p) = E sym (p ) 



(5) 



where u = p/po is the reduced density and E sym (po) = 35 
MeV is the symmetry energy at normal nuclear matter 
density. The symmetry potential acting on a nucleon 
derived from the above nuclear symmetry energy is [24[ 

V sym (p,S) = ±2[£ sym (p V - 12.7u 2 / 3 }5 

+ [£ sym (p )( 7 - l)u 7 + 4.2« 2 / 3 ]«S 2 , (6) 

where "+" and "-" are for neutrons and protons, respec- 
tively. 

For the isospin-independent part of the nuclear EOS, 
E{p, 8 = 0), we use a Skyrme-like parameterization, i.e., 

E{p,5 = Q) = \eW* + ~u+-^-u°, (7) 
5 2 l + o" 

where the first term is the kinetic part and E F is the 
Fermi energy of symmetric nuclear matter at normal nu- 
clear matter density. The parameters a = —358.1 MeV, 
b = 304.8 MeV, and a = 7/6 correspond to the so-called 
soft EOS with incompressibility Kq = 201 MeV, while 
a = -123.6 MeV, b = 70.4 MeV, and a = 2 give the so- 
called stiff EOS with incompressibility K = 380 MeV. 
Both EOS's have a saturation density po = 0.16 fin -3 
for symmetric nuclear matter. The isospin-independent 
potential obtained from this EOS is then 
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The pressure of nuclear matter is given by 

2 dE(p,S) 



P{p,S) 



dp 



(9) 



which leads to a contribution from nuclear symmetry en- 
ergy given by 

P sym (p,<5) = ^ ggp(g) 

= PQ E sym (pohS 2 u'< + \ (10) 
besides that due to isospin-independent energy given by 

o ba 



Po(p,S = 0) = Ie° fPo u 5 / 3 + 1 



apou 



1 + a 



(11) 



At moderate and high temperatures T (T > 4 MeV), 
the pressure can also be decomposed into an isospin- 
dependent part P sym (p, 6, T) and an isospin-independent 
part P {p,S = 0,T), i.e., 



oo / \ 3 \ n 



+ (l-S) 1+n -2}+ P0 S 2 

^ S ym(po)7« 7+1 " 8.5 U 5 / 3 1 , (12) 



and 



Po{p,S 
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o,r) = Tp[i + ^6 n (^r] 
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where &^s are the inversion coefficients given in Ref. 0, 
and 
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is the thermal wavelength of nucleons with mass m. In 
Eq. (12), we have used the explicit ^-dependent kinetic 
part of symmetry energy instead of the parabolic approx- 
imation to calculate the pressure at finite temperatures. 

In the present work, we consider two cases of 7 = 0.5 
(soft) and 2 (stiff) to explore the large range of L and 
K sym discussed above. The 7 = 0.5 (2) gives L = 52.5 
(210.0) MeV and K sym = -78.8 (630.0) MeV. In Figs. □ 
(a) and (b), we show the nuclear symmetry potentials for 
protons and neutrons in isospin asymmetric nuclear mat- 
ter with S = 0.2 and 0.4, respectively. Except at densities 
below 0.5po for the case of 7 = 2, the symmetry poten- 
tial is attractive for protons and repulsive for neutrons 
and increases with density and isospin asymmetry. Fur- 
thermore, the stiff symmetry potential becomes stronger 
than the soft one when nuclear density exceeds certain 
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FIG. 1: (Color online) Neutron and proton symmetry poten- 
tials as functions of density in asymmetric nuclear matter with 
isospin asymmetry (a) <5 = 0.2 and (b) S = 0.4 for both stiff 
(7 = 2) and soft (7 = 0.5) symmetry energies, respectively. 
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FIG. 2: (Color online) Pressure of asymmetric nuclear matter 
due to the symmetry energy, P sy m, as a function of density 
for different temperatures and isospin asymmetries: (a) T — 
MeV and 5 = 0.2, (b) T = 10 MeV and S = 0.2 (c) T = 
MeV and 8 = 0.4, and (d) T = 10 MeV and S = 0.4. Solid 
and dashed curves are for soft (7 = 2.0) and stiff (7 = 0.5) 
symmetry energies, respectively. 



value around normal nuclear density, that depends on 
the isospin asymmetry and is different for neutrons and 
protons. 

Different symmetry energies also result in different 
thermodynamical properties for neutron-rich matter. 
The pressure due to symmetry energy, P sym , in an asym- 
metric nuclear matter with isospin asymmetry d = 0.2 or 
0.4 and temperature T = or 10 MeV, is shown in Fig. 
121 as a function of density for both soft (7 = 0.5) and 
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FIG. 3: (Color online) Isospin-independent part of pressure, 
Po, as a function of density for T — (a) and 10 (b) MeV, 
respectively, with Ko= 380 (solid curves) and 201 (dashed 
curves) MeV. 



stiff (7 = 2) symmetry energies. It is seen that the stiff 
symmetry energy gives a larger pressure in asymmetric 
nuclear matter than the soft one at same density, and 
the difference increases with both density and isospin 
asymmetry. Since the symmetry energy is insensitive 
to temperature, temperature dependence of the pressure 
from symmetry energy is rather weak as in more sophis- 
ticated finite-temperature Skyrme-Hartree-Fock calcula- 
tions |49j. For comparison, the isospin-independent part 
of pressure, Po, is shown in Fig. [31 as a function of den- 
sity for T = (a) and 10 (b) MeV, and with K = 201 
(dashed curves) and 380 (solid curves) MeV. We see that 
the stiff EOS gives a stronger pressure at both higher 
densities and temperatures. 

Although the total incompressibility of asymmetric nu- 
clear matter is mainly determined by the stiffness of the 
isoscalar part of nuclear energy, the stiffness of symmetry 
energy has non-negligible effects on its EOS and thermo- 
dynamical properties. We thus expect that different sym- 
metry energies would lead to different space-time evo- 
lutions in heavy-ion collisions induced by neutron-rich 
nuclei at intermediate energies, resulting in observable 
effects in the nucleon-nucleon correlation functions. 



III. ISOSPIN-DEPENDENT 
BOLTZMANN-UEHLING-UHLENBECK (IBUU) 
TRANSPORT MODEL 



Our study is based on an isospin-dependent 
Boltzmann-Uehling-Uhlenbeck (IBUU) transport 
model (e.g., [13, 12a, |2fj |5(J, |51| ) In this model, initial 
positions of protons and neutrons are determined accord- 
ing to their density distributions from the relativistic 
mean-field (RMF) theory HHHQ. The neutron and 
proton initial momenta are then taken to have uniform 
distributions inside their respective Fermi spheres with 



Fermi momenta determined by their local densities via 
the Thomas-Fermi approximation. For the isoscalar 
potential, we use as default the Skyrme potential 
with incompressibility Ko — 380 MeV. Although the 
transverse flow data from heavy ion collisions are best 
described by a momentum-dependent soft potential with 
A'o = 210 MeV [oil l56|. they can also be approximately 
reproduced with a momentum-independent stiff poten- 
tial with Ko = 380 MeV. Dependence of the reaction 
dynamics on isospins is included through the isospin- 
dependent total and differential nucleon-nucleon cross 
sections and Pauli blockings, the symmetry potential 
V synu and the Coulomb potential V c for protons. For 
nucleon-nucleon cross sections, we use as default the 



experimental values in free space <j e 



For a review 



of the IBUU model, we refer the reader to Ref. |3j. 
To study effects due to the medium dependence of 
nucleon-nucleon cross sections, we also use the pa- 
rameterized in-medium nucleon-nucleon cross sections 
fin-medium fr° m the Dirac-Brueckner approach based on 
the Bonn A potential [57j |. In Figs. E]( a ), (b) and (c), we 
show the density dependence of total nucleon-nucleon 
cross sections cr exp (squares) and cri n - m edium (lines) for 
laboratory incident energies E = 50, 100 and 200 MeV, 
respectively. In experimental free-space cross sections, 
the neutron-proton cross section is about a factor of 3 
larger than the neutron-neutron or proton-proton cross 
section. On the other hand, in-medium nucleon-nucleon 
cross sections from the Dirac-Brueckner approach based 
on the Bonn A potential have smaller magnitude and 
weaker isospin dependence than tr exp but strong density 
dependence. 
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FIG. 4: (Color online) Density dependence of total nucleon- 
nucleon cross sections at laboratory energies E — 50 (a), 100 
(b) and 200 (c) MeV, for the experimental free-space cross sec- 
tions (squares) and the in-medium nucleon-nucleon cross sec- 
tions from the Dirac-Brueckner approach based on the Bonn 
A potential (lines). 

To solve the IBUU model, we use the usual test par- 
ticle method |5fJ. Results presented in the following are 
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obtained with 10, 000 events using 100 test particles for 
a physical nucleon. 



IV. NUCLEON EMISSION FUNCTIONS 

The nucleon emission function describes the distribu- 
tions of time, position, and momentum at which nucleons 
are emitted. In the following, we first show the distribu- 
tion of nucleon emission times and the dependence of the 
average emission time on nucleon momentum. This is 
then followed by the distribution of nucleon emission po- 
sitions and its dependence on nucleon momentum. The 
momentum distribution of emitted nucleons is also given. 
All momenta are in the center-of-mass system unless 
stated otherwise. 



A. Nucleon emission times 

The emission times of different particles in heavy ion 
collisions are relevant for understanding both the colli- 
sion dynamics and the mechanism for particle produc- 
tion. In heavy ion collisions at intermediate energies, 
the dynamics of nucleon emissions is mainly governed 
by the pressure of the excited nuclear matter produced 
during the initial stage of collisions 0> EH • Since the 
stiff symmetry energy gives a larger pressure than that 
due to the soft symmetry energy as shown in Fig. [21 it 
leads to a faster emission of neutrons and protons. The 
relative emission sequence of neutrons and protons is, 
however, determined by the difference in their symmetry 
potentials. With the symmetry potential being generally 
repulsive for neutrons and attractive for protons, neu- 
trons are expected to be emitted earlier than protons. 
The difference in neutron and proton emission times is 
particularly large for the soft symmetry potential as its 
magnitude at densities below the normal nuclear matter 
density, where most nucleons are emitted, is larger than 
that of stiff symmetry potential as seen in Fig. ^ Since 
the stiff symmetry potential changes sign when nuclear 
density drops below 0.5po, the relative emission sequence 
of neutrons and protons in this case depends on the aver- 
age density at which they are emitted. Also, the emission 
time for protons is affected by the repulsive Coulomb po- 
tential. Of course, details of these competing effects on 
nucleon emissions depend on both the reaction dynam- 
ics and momenta of nucleons. These details can only be 
studied by using transport models. 

As an example, we study here central collisions of 52 Ca 
+ 48 Ca at E = 80 MeV/nucleon. This particular reac- 
tion system with isospin asymmetry 5 — 0.2 can be stud- 
ied at RIA. In the present work, nucleons are considered 
as being emitted when their local densities are less than 
po/8 and subsequent interactions does not cause their 
recapture into regions of higher density. Other emission 
criteria, such as taking the nucleon emission time as its 
last collision time in the IBUU model, do not change 
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FIG. 5: (Color online) Emission rates of protons and neu- 
trons as functions of time for different cases: (a) Soft or stiff 
symmetry energy with Ko= 380 MeV and a cxp ; (b) different 
EOS's, i.e., 7V = 380 and 200 MeV with 7 = 0.5; and (c) 
different JV-JV cross sections with Kq= 380 MeV and 7 = 0.5. 



our conclusions. In Fig. we show the emission rates 
of protons and neutrons as functions of time for differ- 
ent cases. It is clearly seen that there are two stages of 
nucleon emissions; an early fast emission and a subse- 
quent slow emission. This is consistent with the long- 
lived nucleon emission source observed in previous BUU 
calculations [53. In Fig. [S] (a), effects due to differ- 
ent symmetry energies are shown, and we find that the 
stiff symmetry energy (triangles) enhances the emission 
of early high momentum protons and neutrons but sup- 
presses late slow emission compared with results from 
the soft symmetry energy (squares). Difference between 
the emission rates of protons and neutrons is, however, 
larger for the soft symmetry energy. Both features are 
what we have expected from above discussions about the 
effects of nuclear symmetry energy on particle emissions. 
Fig. [S] (b) shows results by using different incompress- 
ibilitics of Kq — 380 (squares) and 201 (triangles) MeV 
for the isoscalar potential. It is seen that the nucleon 
emission rate is not sensitive to the incompressibility Kq 
of the symmetric nuclear matter EOS, except that the 
stiff EOS slightly enhances the nucleon emission rate at 
early stage of reactions. In Fig. 0(c), we compare results 
from the experimental free space JV-JV cross section <7 exp 
(squares) and those from the in-medium JV-JV cross sec- 
tions <7i n _ me dium (triangles). Again, the nucleon emission 
rate is found to be insensitive to the JV-JV cross sections 
used for the present reaction system, except that the in- 
medium cross sections slightly enhance the nucleon emis- 
sion rate at later stage of reactions (50 ~ 100 fm/c). 
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FIG. 6: Time evolutions of scaled quadrupole moment of nu- 
cleoli momentum distribution Q zz /A Tes by using Ko = 380 
and 210 MeV with 7 = 0.5, Ko = 380 with 7 = 2, and 
Ko — 380 MeV with 7 = 0.5 and in-medium JV-JV cross 
sections. 
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FIG. 7: Time evolution of average central density by using 
K = 380 and 210 MeV with 7 = 0.5, K = 380 with 7 = 2, 
and Ko = 380 MeV with 7 = 0.5 and in-medium JV-JV cross 
sections. 



To see if nucleons are emitted from an equilibrium or 
non-equilibrium sources, we consider time evolution of 
the scaled quadrupole moment of nucleon momentum dis- 
tribution Q zz /A rcB denned as [60| 

Q zz /A ICS (t) = l/A res J ^L(2pl~pl-p 2 y )f(r,p,t), 

(15) 

where A rcs is the mass number of the residue composed 
of nucleons at local density larger than p /8 and /(r, p, t) 
is the nucleon phase-space distribution function given by 
the IBUU model. How far the emission source devi- 
ates from thermal equilibrium is then given by the value 
of Q zz , with Q zz — implying that thermal equilib- 
rium is achieved. Shown in Fig. [5] are time evolutions 
of Qzz/A res for different nuclear symmetry potentials, 
isoscalar potentials, and JV-JV cross sections. It is seen 
that the emission source nearly reaches thermal equi- 
librium after around 40 fm/c except the case with in- 
medium JV-JV cross sections (dash-dotted curve). In the 
latter case, the thermal equilibration time is delayed to 
about 70 fm/c as a result of smaller in-medium N-N cross 
sections than free N-N cross sections. Therefore, early 
emitted nucleons are mostly from pre-equilibrium stage 
of heavy ion collisions, while late emitted ones are more 
like statistical emissions. Since positive Q zz indicates in- 
complete nuclear stopping or nuclear transparency while 
negative one implies transverse expansion or collectivity, 
all cases except the one with in-medium N-N cross sec- 
tions thus show strong nuclear stopping but weak trans- 
verse expansion after 120 fm/c. 

The weak dependence of the nucleon emission rate on 
the stiffness of isoscalar potential is due to reduced dif- 
ference in the pressure of excited nuclear matter as a re- 
sult of different maximum densities reached in collisions, 
with the stiff one giving a lower density than the soft 



one. This can be seen in Fig. [7] where the time evo- 
lution of average central density (calculated in a sphere 
located at the center of colliding nuclear matter with ra- 
dius of 3 fm) is shown by using K = 380 (solid curve) or 
201 (dashed curve) MeV but with same soft symmetry 
energy and free nucleon- nucleon cross sections. Maxi- 
mum average central densities reached in the collisions 
are about 1.46po and 1.66po, respectively, for K = 380 
and 201 MeV, resulting in similar pressures of about 8 
and 7 MeV/fm 3 , respectively. In addition, the soft EOS 
leads to a slower expansion and a longer higher density 
stage than the stiff EOS. On the other hand, the maxi- 
mum average central density reached in the collisions is 
reduced only slightly if the stiff symmetry energy is used 
for the case of Ko = 380 MeV and free JV-JVcross sections 
(dotted curve). Also shown in Fig. 0is the time evolu- 
tion of average central density by using Ko — 380 MeV 
and the soft symmetry energy but with in-medium N-N 
cross sections <7i n _ me di U m (dash-dotted curve). One finds 
that different JV-JV cross sections lead to almost similar 
time evolutions of average central density. Although the 
free JV-JV cross sections, which have larger values than 
the in-medium ones, would lead to more stopping and 
thus higher density, enhanced Pauli blocking and mean- 
field pressure effects bring the resulting maximum density 
similar to that obtained using the in-medium cross sec- 
tions. As a result, the two different cross sections lead to 
similar pressures during the high density stage of heavy 
ion collisions, and thus give comparable nucleon emission 
times as shown in Fig. [S] (c) . 

Energies of emitted nucleons are correlated to the times 
at which they are emitted. Generally, earlier emitted nu- 
cleons have higher energies than later emitted ones. This 
can be seen in Fig. [S] where we show the average emis- 
sion times of protons and neutrons as functions of their 
momenta by using Ko = 380 MeV with the soft or stiff 
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FIG. 8: (Color online) Average emission times of protons and 
neutrons as functions of their momenta for different symmetry 
energies. 



symmetry energy. Indeed, high momentum nucleons are 
emitted in the early pre-equilibrium stage of collisions 
when the average density is relatively high, while low 
momentum ones are mainly emitted when the system is 
close to equilibrium and the average density is low. Fig. 
[S] further shows that the average emission time of nucle- 
ons with a given momentum is earlier for the stiff sym- 
metry energy (triangles) than for the soft one (squares). 
Moreover, there are significant delays in proton emissions 
(filled squares) with the soft symmetry energy. These fea- 
tures are consistent with the results shown in Fig. OJa). 
For stiff symmetry energy, the change in the emission se- 
quence of neutrons and protons with their momenta is 
consistent with the change in the sign of stiff symmetry 
potential at low densities as shown in Fig. 

We note, however, the distribution of the emission 
times for nucleons with same momentum is quite broad 
particularly for low momentum nucleons. This is shown 
in Fig. |5|for nucleons with momenta of 100 (a), 200 (b), 
and 300 (c) MeV/c (all with bin size of 20 MeV/c), ob- 
tained with Kq = 380 MeV, free JV-JV cross sections, and 
either soft or stiff symmetry energy. Although the peak 
emission time increases with decreasing nucleon momen- 
tum, the distribution of nucleon emission times is broad. 
Furthermore, at later times more nucleons with momen- 
tum of 200 MeV/c are emitted than those with momen- 
tum of 100 MeV/c are emitted. For both protons and 
neutrons, more are emitted for the stiff than for the soft 
symmetry energy. 

In low and intermediate energy heavy-ion collisions in- 
duced by heavy nuclei, the Coulomb interaction is ex- 
pected to affect the emission time of protons and thus 
also the neutron-proton correlation function |6lL l62j | . In 
present IBUU simulations, the Coulomb potential is in- 
cluded explicitly in the dynamic evolution of colliding 
nuclei. Effects of Coulomb interaction on the proton 
emission time can be seen in Fig. 1101 where the aver- 
age emission time of protons is shown as a function of 
their momenta with (filled squares) and without (open 
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FIG. 9: (Color online) Emission rates of neutrons and protons 
with momenta 100 (a), 200 (b), and 300 (c) MeV/c (with bin 
size of 20 MeV/c) for Ko = 380 MeV, free JV-JV cross sections, 
and either soft or stiff symmetry energy. 
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FIG. 10: (Color online) Average emission time of protons as 
a function of their momenta with (filled squares) and without 
(open squares) Coulomb interaction by using the soft (a) or 
stiff (b) symmetry energy. 



squares) Coulomb interaction by using either the soft 
(left panel) or stiff (right panel) symmetry energy. One 
finds that the Coulomb interaction is not important for 
52 Ca -I- 48 Ca as it only shortens the proton average emis- 
sion time by factors less than 10% for both the soft and 
stiff symmetry energies. The Coulomb potential effect 
is even smaller for protons with high momentum, which 
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are emitted early during collisions. We note that the 
Coulomb potential does not affect the neutron emission 
time as expected. 
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FIG. 11: (Color online) Root-mean-squared transverse emis- 
sion radius (ii rms ) and longitudinal emission radius (Z ims ) of 
neutrons and protons as functions of their momenta by using 
Ko ~ 380 MeV with either soft (squares) or stiff (triangles) 
symmetry energy. 

The size of emission source can be characterized by the 
root-mean-squared radii of emitted protons and neutrons 
in the transverse ((i? r ms)) and longitudinal ((Z rms )) di- 
rections. In Fig. ^2 we show these radii as functions of 
nucleon momentum using Ko = 380 MeV and free N-N 
cross sections but with either soft (squares) or stiff (trian- 
gles) symmetry energy. It is seen that for both symmetry 
energies the emission source has an oblate shape with a 
larger transverse radius than longitudinal radius. Fur- 
thermore, the emission source size is larger for the soft 
than the stiff symmetry energy. However, the difference 
between the transverse radii of proton and neutron emis- 
sion sources is more appreciable in the case of soft sym- 
metry energy than stiff symmetry energy. These features 
are consistent with the emission times shown in Fig. [5] 

The decrease of the transverse radius with increasing 
nucleon momentum shown in Fig. 1111 is expected as 
higher momentum nucleons are emitted earlier when the 
source size is more compact. The small transverse radius 
of the emission source for lower momentum nucleons are, 
on the other hand, due to the broad distributions in nu- 
cleon emission times, resulting in the emissions of more 
intermediate momentum nucleons than low momentum 
nucleons at later times. As a result, the source size for 
low momentum nucleons is largely determined by ear- 
lier emitted nucleons. This can also been seen in Fig. 
1121 where we show the transverse radius distributions of 
emitted nucleons with momenta 100 (a), 200 (b), and 
300 (c) MeV/c (with bin size of 20 MeV/c) for K = 380 



FIG. 12: (Color online) Transverse radius distribution of 
emitted nucleons with momenta 100 (a), 200 (b), and 300 (c) 
MeV/c (with bin size of 20 MeV/c) for K = 380 MeV, free 
N-N cross sections, and either soft or stiff symmetry energy. 



MeV, free N-N cross sections, and either soft or stiff 
symmetry energy. It is seen that the peak of the trans- 
verse radius distribution first increases then decreases as 
the nucleon momentum gets larger, consistent with that 
shown in Fig. ^] Furthermore, width of the radius dis- 
tribution is narrower for high momentum nucleons than 
for low momentum ones. 
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FIG. 13: Density contours p(x, 0, z) at different times by 
using Ko = 380 MeV and free N-N cross sections with the 
soft symmetry energy. 

To see more clearly the size and structure of nucleon 
emission source, we show in Figs. 1131 and 1141 the den- 
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sity reached in collisions. In Fig. ^] we show the time 
evolution of maximum density for different incompress- 
ibilitics, N-N cross sections, and symmetry energies as in 
Fig. [7| We see that after about 30 fm/c when the bubble 
is formed the maximum density has a value around nor- 
mal nuclear density and is very different from the central 
density shown in Fig. Since low momentum nucleons 
are mainly emitted during the late stage of heavy ion 
collisions from the ring and its decaying fragments, the 
symmetry energy effect is thus small as a result of the 
small difference between soft and stiff symmetry poten- 
tials at normal nuclear density as shown in Fig. 



FIG. 14: (Color online) Density contours p(x, y, 0) at dif- 
ferent times by using Ko = 380 MeV and free N-N cross 
sections with the soft symmetry energy. The thick curves 
represent po/8. 
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FIG. 15: Time evolution of maximum density for different in- 
compressibilities, N-N cross sections, and symmetry energies 
as in Fig. Q 



C. Momentum distributions of emitted nucleons 

Fig. 1161 shows the momentum distributions of emitted 
neutrons and protons for K — 380 MeV and free N-N 
cross sections with soft (squares) or stiff (triangles) sym- 
metry energy. It is seen that most nucleons are emitted 
with momenta around p = 180 MeV/c, and the peak mo- 
mentum shifts to a larger value as the symmetry energy 
becomes softer. The lower peak momentum in the case 
of stiff symmetry energy is due to the larger potential en- 
ergy during initial compressions and thus lower average 
nuclcon kinetic energy as a result of energy conservation. 
Furthermore, symmetry energy effects on lower momen- 
tum nucleons are stronger for protons than for neutrons, 
which is consistent with the larger variation of proton 
symmetry potential at low densities with the stiffness of 
symmetry energy than that of neutron symmetry poten- 
tial as shown in Fig. 



sity contours p(x, 0, z) in the reaction plane and p{x, y, 0) 
perpendicular to the reaction plane at different times by 
using Kq = 380 MeV and free N-N cross sections with 
the soft symmetry energy. It is seen that after initial con- 
tact of target and project at 1 fm/c, the nuclear matter 
is compressed to high density during the first 20 fm/c. 
This is followed by expansions, mainly in the transverse 
direction, leading to the formation of an oblate bub- 
ble nuclear matter, which subsequently changes into a 
ring-shaped structure. The ring structure has a long 
lifetime and decays into several fragments after about 
100 fm/c. Such exotic structures also exist in the other 
cases, i.e., using Ko = 201 MeV, stiff symmetry energy, 
or in-medium N-N cross sections. This interesting phe- 
nomenon is not new and was studied extensively a decade 
ago. For medium mass nuclei, this is a common phe- 
nomenon in central or near central collisions at energy we 
consider here |U l6"H l6l| . The formation of bubble and 
ring structures during expansion implies that the central 
density shown in Fig. [7] is not always the maximum den- 




(MeV/c) 



FIG. 16: (Color online) Momentum distributions of emitted 
neutrons and protons for Ko = 380 MeV and free N-N cross 
sections with soft (squares) or stiff (triangles) symmetry en- 
ergy. 
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V. NUCLEON-NUCLEON CORRELATION 
FUNCTIONS 

The nudeon emission function, which is important for 
understanding the reaction dynamics in heavy ion col- 
lisions, can be extracted from two-particle correlation 
functions; see, e.g., Refs. [27l l28l 123. l30| for reviews. 
In most studies, only the two-proton correlation function 
has been measured [5^, |6^, ||3]l EH E3 • Recently, data 
on two-neutron and neutron-proton correlation functions 
have also become available. The neutron-proton correla- 
tion function is especially useful as it is free of correlations 
due to wave-function anti-symmetrization and Coulomb 
interactions. Indeed, Ghctti et al. have deduced from 
measured neutron-proton correlation function the time 
sequence of neutron and proton emissio ns FpLFj) . 

In standard Koonin-Pratt formalism [72 

flzllzi, the 

two-particle correlation function is obtained by convolut- 
ing the emission function g(jp,x), i.e., the probability for 
emitting a particle with momentum p from the space- 
time point x = (r,i), with the relative wave function of 
the two particles, i.e., 

r ,p s _ Jd 4 x 1 d i x 2 g(P/2,x 1 )g(-p/2,x 2 )\ ( j ) (q,T)\ 2 
[ ,qj f d*x ig (P/2,xi) f d*x 2 g(P/2,x 2 ) ■ 

(16) 

In the above, P(= pi + P2) and q(=|(pi — P2)) are, re- 
spectively, the total and relative momenta of the particle 
pair; and 0(q, r) is the relative two-particle wave function 
with r being their relative position, i.e., r = (r 2 — ri) — 
|(vi + v 2 )(^2 — t\)- This approach has been very use- 
ful in studying effects of nuclear equation of state and 
nucleon-nucleon cross sections on the reaction dynam- 
ics of intermediate energy heavy-ion collisions J28j. In 
the present paper, we use the Koonin-Pratt method to 
determine the nucleon-nucleon correlation functions in 
order to study the effect due to the density dependence 
of nuclear symmetry energy on the spatial and tempo- 
ral structure of nucleon emission source in intermediate 
energy heavy ion collisions. 

A. Symmetry energy effects 

Using the program Correlation After Burner [75j . 
which takes into account final-state nucleon-nucleon in- 
teractions, we have evaluated two-nucleon correlation 
functions from the emission function given by the IBUU 
model. Shown in Fig. 1171 are two-nucleon correla- 
tion functions gated on total momentum P of nucleon 
pairs from central collisions of 52 Ca + 48 Ca at E = 80 
MeV/nucleon. The left and right panels are for P < 
300 MeV/c and P > 500 MeV/c, respectively For 
both neutron-neutron (upper panels) and neutron-proton 
(lower panels) correlation functions, they peak at q rs 
MeV/c. The proton-proton correlation function (mid- 
dle panel) is, however, peaked at about q = 20 MeV/c 
due to strong final-state s-wave attraction, but is sup- 



pressed at q = as a result of Coulomb repulsion and 
anti-symmetrization of two-proton wave function. These 
general features are consistent with those observed in ex- 
perimental data from heavy ion collisions [7C|. 
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FIG. 17: (Color online) Two-nucleon correlation functions 
gated on total momentum of nucleon pairs using the soft 
(filled squares) or stiff (filled triangles) symmetry energy. Left 
panels are for P < 300 MeV/c while right panels are for 
P > 500 MeV/c. 

Since emission times of low- momentum nucleons do not 
change much with the different E sym (p) used in IBUU 
model as shown in Fig. |SJ two-nucleon correlation func- 
tions are not much affected by the stiffness of symme- 
try energy. On the other hand, the emission times of 
high-momentum nucleons, which are dominated by those 
with momenta near 250 MeV/c, differ appreciably for 
the two symmetry energies considered here. Correlation 
functions of high-momentum nucleon pairs thus show a 
strong dependence on nuclear symmetry energy. Gat- 
ing on nucleon pairs with high total momentum allows 
one to select those nucleons that have short average spa- 
tial separations at emissions and thus exhibit enhanced 
correlations. For these nucleon pairs with high total mo- 
mentum, strength of their correlation function is stronger 
for the stiff symmetry energy than for the soft symmetry 
energy: about 30% and 20% for neutron-proton pairs and 
neutron-neutron pairs at low relative momentum q = 5 
MeV/c, respectively, and 20% for proton-proton pairs at 
q = 20 MeV/c. The neutron-proton correlation func- 
tion thus exhibits the highest sensitivity to variations of 
E sym {p)- As shown in Fig. |S|and discussed earlier, the 
emission sequence of neutrons and protons is sensitive 
to E sym (p). With stiff E sym (p), high momentum neu- 
trons and protons are emitted almost simultaneously, and 
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they are thus temporally strongly correlated, leading to a 
larger neutron-proton correlation function. On the other 
hand, proton emissions in the case of soft E sym (p) are 
delayed compared to neutrons, so they are temporally 
weakly correlated with neutrons. The resulting values 
of neutron-proton correlation function are thus smaller. 
Furthermore, both neutrons and protons are emitted ear- 
lier with stiff E sym (p), so they have shorter spatial sep- 
aration at emissions. As a result, neutrons and protons 
are more correlated for the stiff symmetry energy than 
for the soft symmetry energy. Our results thus clearly 
demonstrate that correlation functions of nucleon pairs 
with high total momentum can indeed reveal sensitively 
the effect of nuclear symmetry energy on the temporal 
and spatial distributions of emitted nucleons. 



ever, the effect becomes weaker when their total momen- 
tum becomes very large. This is due to the fact that very 
energetic neutrons and protons mainly come from very 
early stage of collisions when the reaction system reaches 
maximum compression and nucleon emissions are largely 
affected by the isospin-independent part of EOS. The 
symmetry energy effect on the emission time of nucle- 
ons is therefore reduced for these very energetic neutrons 
and protons as shown in Fig. [SJ Furthermore, emission 
times of very energetic neutrons and protons are strongly 
influenced by the high density behavior of symmetry en- 
ergy, and for the stiff symmetry energy very high momen- 
tum neutrons are emitted earlier than protons as shown 
in Fig. [S] As a result, the symmetry energy effect on 
the correlation function of very energetic neutron-proton 
pairs is also reduced. 
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FIG. 18: (Color online) Total momentum dependence of the 
correlation functions for neutron-neutron (a) and neutron- 
proton (c) pairs with relative momentum q — 5 MeV/c and 
proton-proton (b) pairs with relative momentum q = 20 
MeV/c by using Ko = 380 MeV with both soft (filled squares) 
and stiff (open squares) symmetry energies. 

To see more clearly the symmetry energy effect on 
nucleon-nucleon correlation functions, we show in Fig. 
1181 the total momentum dependence of the correlation 
functions for neutron-neutron and neutron-proton pairs 
with relative momentum q — 5 MeV/c and proton- proton 
pairs with relative momentum q = 20 MeV/c by using 
Kq = 380 MeV with both soft (filled squares) and stiff 
(open squares) symmetry energies. It is seen that the 
symmetry energy effect is more pronounced for pairs with 
higher total momentum as discussed previously. How- 
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FIG. 19: (Color online) Same as Fig. 1 171 but with angular cut 
0c.m. = 90° ± 10° in center-of-mass system. 

Since the nuclear symmetry energy affects differently 
the transverse and longitudinal radii of nucleon emission 
source, it is useful to study also the dependence of two- 
nucleon correlation functions on the direction of emit- 
ted nucleons. We have thus considered both transverse 
cut with 9 c m — 90° ± 10° and longitudinal cut with 
c . m . = 0° ~ 20° or 160° - 180°, where c . m . is the 
center-of-mass angle with respect to the beam direction. 
Results similar to Fig. 1171 for two-proton, two-neutron, 
and neutron-proton correlation functions using both stiff 
and soft symmetry energies but stiff isospin independent 
EOS and free JV-JV cross sections are shown in Fig. IT^l 
and Fig. [5D|for transverse and longitudinal cuts, respec- 
tively. It is seen that the transverse cut gives smaller 
values for nucleon-nucleon correlation functions than the 
longitudinal cut. This is easy to understand as the corre- 
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FIG. 20: (Color online) Same as Fig. [13 but w ith angular cut t (fm/c) 

9cm.. = 0° ~ 20° or 160° ~ 180° in center-of-mass system. 



lation functions with transverse (longitudinal) cut mainly 
reflect the transverse (longitudinal) size or lifetime of 
emission source. A transverse radius larger than the lon- 
gitudinal radius shown in Fig. ^]thus leads to smaller 
values for the nucleon-nucleon correlation functions. Our 
results further show that the correlation functions of nu- 
cleon pairs with total momentum P > 500 MeV/c ob- 
tained with transverse cut have stronger symmetry en- 
ergy dependence than those obtained with longitudinal 
cut. For transverse cut, the symmetry energy effects are 
about 23% and 33% for neutron-neutron and neutron- 
proton pairs at q — 5 MeV/c, and 22% for proton- proton 
pairs at q = 20 MeV/c, respectively, but corresponding 
effects are 14%, 22%, and 19%, respectively, for longitu- 
dinal cut. However, compared to results in Fig. I17l with- 
out the angular cut, there is only a slight enhancement 
of two-nucleon correlation functions when the transverse 
cut is imposed, as most nucleons are emitted in the trans- 
verse direction. 

In the above analysis, the IBUU calculations are ter- 
minated at 200 fm/c. In Fig. [21 we show the stop time 
dependence of the correlation functions for high total mo- 
mentum (P > 500 MeV) neutron-neutron and neutron- 
proton pairs with relative momentum q = 5 MeV/c and 
proton-proton pairs with relative momentum q = 20 
MeV/c by using K = 380 MeV with soft (filled squares) 
or stiff (open squares) symmetry energy. It is seen that 
values of these correlation functions are essentially fixed 
after about 150 fm/c for both soft (open squares) and 
stiff (filled squares) symmetry energies. Furthermore, the 
symmetry energy effect starts to appear after 50 fm/c, 
when the central density is about l/5po of normal nuclear 



FIG. 21: (Color online) Dependence of two-nucleon cor- 
relation functions on stop times for high total momentum 
(P > 500 MeV) neutron-neutron (a) and neutron-proton (c) 
pairs with relative momentum q — 5 MeV/c and proton- 
proton (b) pairs with relative momentum q = 20 MeV/c by 
using Kq = 380 MeV with soft (filled squares) or stiff (open 
squares) symmetry energy. 

matter density as shown in Fig. and the maximum 
density is about l/2po as shown in Fig. 1151 indicating 
that two-nucleon correlation functions are sensitive to 
the properties of nuclear symmetry energy at densities 
around and below the normal nuclear density. 



B. EOS and N-N cross section effects 

We have also studied the effects of isoscalar potential 
and in-medium N-N cross sections on nucleon-nucleon 
correlation functions. In Fig. [22 we show gated two- 
nucleon correlation functions from central collisions of 
52 Ca + 48 Ca at B = 80 MeV/nucleon using the soft 
symmetry energy but different values of incompressibil- 
ity Kq and N-N cross sections. One finds that reducing 
K from 380 (filled squares) to 201 (filled triangles) MeV 
increases the value of the correlation function for high 
momentum pairs (P > 500 MeV/c) but decreases that 
for low momentum pairs (P < 300 MeV/c). Specifically, 
the increases are about 8% and 13% for high momen- 
tum neutron-proton and neutron-neutron pairs with rel- 
ative momentum q — 5 MeV/c, and about 5% for high 
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FIG. 22: (Color online) Two-nucleon correlation functions 
gated on total momentum of nucleon pairs using the soft sym- 
metry energy but different values of incompressibility Ko and 
JV-JV cross sections. Left panels are for P < 300 MeV/c while 
right panels are for P > 500 MeV/c. 



momentum proton-proton pairs with relative momentum 
q — 20 MeV/c, while corresponding decreases for low mo- 
mentum pairs are about 14%, 3%, and 12%. Changing 
experimental free space JV-JV cross section <7 oxp to in- 
medium JV-JV cross sections Cm-medium (open triangles) 
leads to changes of about 5%, 6%, and 1% in the neutron- 
proton, neutron-neutron, and proton-proton correlation 
functions for high momentum particle pairs (P > 500 
MeV/c), but about 10%, 2%, and 5% for low momentum 
particle pairs (P < 300 Me V/c) . These results are sim- 
ilar to those found in Ref. 28] . The weak dependence 
of two-nucleon correlation functions on the stiffness of 
isoscalar potential and the JV-JV cross sections is consis- 
tent with the nucleon emission times shown in the last 
section. 



C. Impact parameter and incident energy 
dependence 

The nucleon-nucleon correlation functions also depend 
on the impact parameter and incident energy of heavy 
ion collisions. Fig. 1231 shows the impact parameter de- 
pendence of neutron-neutron and neutron-proton corre- 
lation functions at relative momentum q = 5 MeV/c and 
proton-proton correlation function at relative momentum 
q = 20 MeV/c for particle pairs with high total momen- 
tum P > 500 MeV/c by using K = 380 MeV and free 
JV-JV cross sections with the soft (filled squares) or stiff 
(open squares) symmetry energy, ft is seen that the sym- 
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FIG. 23: (Color online) Impact parameter dependence of 
neutron-neutron (a) and neutron-proton (c) correlation func- 
tions at relative momentum q — 5 MeV/c and proton-proton 
(b) correlation function at relative momentum q = 20 MeV/c 
for particle pairs with high total momentum P > 500 MeV/c 
by using A'o = 380 MeV and free JV-JV cross sections with the 
soft (filled squares) or stiff (open squares) symmetry energy. 



metry energy effect is stronger in central and semi-central 
collisions and becomes weaker in semi-peripheral and pe- 
ripheral collisions. This simply reflects the fact that nu- 
clear compression is weaker in peripheral collisions, and 
the symmetry energy effect mainly comes from the small 
difference between the soft and stiff symmetry energies at 
low densities. Also, one finds that values of these corre- 
lation functions decrease with increasing impact param- 
eter, which is consistent with previous BUU results for 
the proton-proton correlation function |6?| . The corre- 
lation function obtained with the stiff symmetry energy 
exhibits, however, a stronger dependence on impact pa- 
rameter than that using the soft symmetry energy as the 
former has a stronger density dependence. 

Dependence of nucleon-nucleon correlation functions 
on the incident energy of heavy ion collisions is shown 
in Fig. |H for high total momentum (P > 500 MeV) 
neutron-neutron and neutron-proton pairs with relative 
momentum q = 5 MeV/c and for proton- proton pairs 
with relative momentum q = 20 MeV/c. The results are 
obtained using Kq = 380 MeV and free JV-JV cross sec- 
tions with the soft (filled squares) or stiff (open squares) 
symmetry energy. It is seen that the values of nucleon- 
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FIG. 24: (Color online) Dependence of nucleon-nucleon corre- 
lation functions on the incident energy for hight total momen- 
tum (P > 500 MeV) neutron-neutron (a) and neutron-proton 
(c) pairs with relative momentum q = 5 MeV/c and proton- 
proton (b) pairs with relative momentum q = 20 MeV/c by 
using Ko = 380 MeV and free N-N cross sections with the 
soft (filled squares) or stiff (open circles) symmetry energy. 



D. Reaction system dependence 
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FIG. 25: (Color online) Two-nucleon correlation functions 
gated on total momentum of nucleon pairs using the soft 
(filled squares) or stiff (filled triangles) symmetry energy 
for the central reaction system 132 Sn + 124 Sn at E= 80 
MeV/nucleon. Left panels are for P < 300 MeV/c while 
right panels are for P > 500 MeV/c. 



nucleon correlation functions increase with increasing in- 
cident energy. This is understandable since nuclear com- 
pression increases with increasing incident energy, and 
nucleons are also emitted earlier, leading to a smaller 
source size. Also, the stiff symmetry energy gives a larger 
value for the nucleon-nucleon correlation function than 
the soft symmetry energy, and the difference does not 
depend much on incident energies, although the relative 
effect is reduced. This is due to the fact that at high in- 
cident energies the symmetry energy effect mainly comes 
from early compression stage of collisions when nuclear 
pressure affects strongly the emission times of nucleons. 
Furthermore, the system expands rapidly at high incident 
energies, the low density behavior of symmetry potential 
thus plays a less important role on the relative emission 
time of neutrons and protons. According to Fig. [8] nu- 
cleons are emitted earlier with stiff rather than with soft 
symmetry energies. Therefore, the symmetry potential 
effects remain appreciable with increasing incident en- 
ergy. On the other hand, the initial compression is rela- 
tively weak at lower incident energies, and the behavior 
of symmetry energy at low densities becomes relevant. 



To see how the above results change for different re- 
action systems, we show in Fig. 1251 two-nucleon correla- 
tion functions for head-on collisions of 132 Sn + 124 Sn at 
E = 80 MeV/nucleon by using the soft (squares) or stiff 
(triangles) symmetry energy. This reaction system has 
a similar isospin asymmetry, i.e., <5 = 0.22 as the reac- 
tion system 52 Ca + 48 Ca studied in the above, and can 
also be studied in future RIA. ft is seen that two-nucleon 
correlation functions in this case are similar to those in 
52 Ca + 48 Ca reactions shown in Fig. El The strength 
of two-nucleon correlation functions at small relative mo- 
menta is, however, about a factor of 2 smaller than that 
in the reaction system 52 Ca + 48 Ca. This reduction is 
mainly due to the larger size of the reaction system 1 Sn 
+ 124 Sn than 52 Ca + 48 Ca. Because of stronger Coulomb 
repulsion for protons, the symmetry energy effect on the 
correlation functions of neutron-proton pairs and proton- 
proton pairs is weaker compared with those in the reac- 
tion system 52 Ca + 48 Ca, i.e., about 17% and 8%, respec- 
tively, instead of 30% and 20%. However, the symmetry 
energy effect on the neutron-neutron correlation function, 
which is not affected by Coulomb potential, is not much 
influenced by the size of reaction system, i.e., it is about 
20% and 15% for the light and heavy reaction systems, 
respectively. As in collisions of 52 Ca + 48 Ca, effects due 
to different isoscalar potentials and JV-JV cross sections 
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on nucleon-nucleon correlation functions are found to be 
small. 



VI. SUMMARY AND OUTLOOK 

In conclusion, using an isospin-dependent transport 
model we explore systematically the isospin effects on 
nucleon-nucleon correlation functions in heavy-ion colli- 
sions induced by neutron-rich nuclei at intermediate en- 
ergies. It is found that the nuclear symmetry energy 
E sy m{p) influences significantly both the emission times 
and the relative emission sequence of neutrons and pro- 
tons. A stiffer E S y m (p) leads to a faster and almost simul- 
taneous emission of neutrons and protons. Consequently, 
two-nucleon correlation functions, especially for neutron- 
proton pairs with high total momentum but low relative 
momentum, are stronger for the stiff symmetry energy 
than the soft one. The correlation function of neutron- 
proton pairs gated on higher total momentum but with 
smaller relative momentum is, on the other hand, insensi- 
tive to the incompressibility of symmetric nuclear matter 
EOS and the in-medium N-N cross sections. This novel 
property can be used as a possible tool to extract use- 
ful information about the density-dependence of nuclear 
symmetry energy. These results do not change much if 
only nucleons emitted in certain directions are consid- 
ered. 

We have also studied the dependence of neutron- 
proton correlation function on the impact parameter and 
incident energy of heavy ion collisions as well as masses of 
the reaction system. We find that the symmetry energy 
effect becomes weaker with increasing impact parame- 
ter and incident energy. Also, the strength of nucleon- 
nucleon correlation function is reduced in collisions of 
heavier reaction systems as a result of larger nucleon 
emission source. For proton-proton and neutron-proton 
correlation functions, the symmetry energy effect is fur- 
ther suppressed by the stronger Coulomb potential in 
heavier reaction systems. Our results thus suggest that 
the correlation function of neutron-proton pairs gated on 



higher total momentum but with smaller relative mo- 
mentum from central or semi-central collisions induced 
by lighter neutron-rich nuclei is a sensitive probe to the 
density dependence of nuclear symmetry energy. We have 
further studied the dependence of symmetry energy effect 
on the time at which IBUU simulations are terminated. 
We find that the symmetry energy effect on nucleon- 
nucleon correlation functions begins to appear when the 
maximum nuclear density is already below the normal 
nuclear matter density. The effect addressed in present 
study via nucleon-nucleon correlation functions is thus 
largely limited to the properties of nuclear symmetry en- 
ergy at subnormal nuclear densities, similar to that from 
studying neutron-skins of radioactive nuclei. 

In the present work, we have not included momen- 
tum dependence in either the isoscalar mean-field po- 
tential or the symmetry potential. The former may af- 
fect the properties of nuclear emission source as it was 
shown that a momentum-dependent mean-field poten- 
tial reduces nuclear stopping or increases nuclear trans- 
parency |7t| . With momentum dependence included in 
the nuclear symmetry potential, differences between the 
neutron and proton potentials are affected not only by 
the density dependence of nuclear symmetry energy but 
also by the magnitude of proton and neutron momenta 
[23. It will be of interest to study how the results ob- 
tained in the present study are modified by the momen- 
tum dependence of nuclear mean field. We plan to carry 
out such a study in the future. 
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